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A HIGH FREQUENCY MODEL OF CASCADE NOISE 


Edmane Envia* 

Acoustics Branch, NASA Lewis Research Center 
Cleveland, Ohio, USA 


ABSTRACT 

Closed form asymptotic expressions for computing high frequency noise generated by 
an annular cascade in an infinite duct containing a uniform flow are presented. There are 
two new elements in this work. First, the annular duct mode representation does not rely on 
the often-used Bessel function expansion resulting in simpler expressions for both the radial 
eigenvalues and eigenfunctions of the duct. In particular, the new representation provides an 
explicit approximate formula for the radial eigenvalues obviating the need for solutions of 
the transcendental annular duct eigenvalue equation Also, the radial eigenfunctions are 
represented in terms of exponentials eliminating the numerical problems associated with 
generating the Bessel functions on a computer. The second new element is the construction 
of an unsteady response model for an annular cascade. The new construction satisfies the 
boundary conditions on both the cascade and duct walls simultaneously adding a new level 
of realism to the noise calculations. Preliminary results which demonstrate the effectiveness 
of the new elements are presented. A discussion of the utility of the asymptotic formulas for 
calculating cascade discrete tone as well as broadband noise is also included. 


INTRODUCTION 

Prediction of fan noise continues to be an 
integral part of the efforts aimed at analysis and 
suppression of turbofan engine noise. To that end, a 
number of schemes have been developed to calculate 
the contribution of various fan-associated noise 
sources. But, whether the goal is to predict rotor inlet 
noise (eg., Mani |1973]), cascade self-noise (eg., 
Glegg [19%]) or rotor-stator interaction noise (e g., 
Envia et. al, [1996]), all of these schemes are 
predicated on computing the acoustic response of a 
blade row to incident flow perturbations 

In strictly non-CFD based methods, which 
continue to be the principal tools of fan noise 
prediction, the blade row acoustic response is 
computed using simplified versions of the cascade 
geometry and/or the governing equations. At the 
simplest end of the spectrum, approximations in both 
the geometry' and governing equations yield analytical 
or semi-analytical expressions for the acoustic 
response. Generally speaking, the approximations 
involve replacing the cascade with zero thickness flat 
plates and linearizing the governing equations about a 
uniform base flow. In this class of methods the 


acoustic response can be determined either directly by 
solving the governing equations or indirectly through 
the use of the annular duct Green’s functions which 
couple the unsteady surface pressure distribution on the 
cascade to the duct acoustic modes. 

However, even at that level of approximation, 
exact solutions are still difficult to obtain, so further 
simplifications are made. Some simplifications are 
based on replacing the annular cascade by a rectilinear 
one at some representative radius and computing the 
acoustic response for the resulting system (eg.. 
Goldstein 1 1976]). In other simplifications, the annular 
geometry is retained but the cascade unsteady pressure 
distribution (henceforth, called the cascade unsteady 
response) is computed at a number of span wise 2D 
strips and the resulting distribution assembled and 
coupled to the duct acoustic modes (e g., Ventres et. al 
[1982]). 

While in both cases useful information 
regarding general trends can be obtained, it is highly 
desirable to retain as much of the cascade and unsteady 
response three-dimensionality as possible Such 3D 
models have been explored by, among others, Namba 
[1977], Kobayashi [1978] and Schulten [1993], but 
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these approaches generally require the use of 
collocation techniques or numerical solutions of partial 
differential equations. 

In this paper, a closed form approximate 
model of cascade noise is presented which takes into 
account the annular geometry of the blade row 
throughout the analysis. The noise field is computed 
by coupling the cascade unsteady response to the duct 
acoustic modes in the standard manner. However, 
asymptotic methods have been employed to circumvent 
the need for numerical solution of the governing 
equations by tailoring the model to the high frequency 
noise regime. Given the current interest in high-speed 
fans, this high frequency specialization, far from being 
a limitation, is in fact an asset since the numerically 
based methods have inherent difficulties handling the 
high frequency response regime 

In what follows, the development of two new 
elements in modeling of annular cascade noise will be 
presented. The first element involves a representation 
of annular duct modes that is not predicated on the use 
of the Bessel functions The second new element 
addresses the construction of an unsteady response to 
incident flow perturbations (i.e., gusts) for an annular 
cascade where the finite-span effects as well as the gust 
three-dimensionality are taken into account. In each 
case, intermediate results will be presented which show 
the effectiveness of the new approximations. The two 
elements are then combined to provide the desired high 
frequency model of cascade noise The paper is 
concluded with a discussion of the utility of the new 
model for predicting cascade tone and broadband 
noise. 


ANALYSIS 

The cascade is modeled as an ensemble of 
zero camber and thickness flat plates enclosed in an 
infinite hard-wall annular duct within which exits a 
uniform and isentropic medium. As shown in Figure 1, 
two coordinate systems will be employed in this work. 
The global system (jc,r,0) is cylindrical polar and is 
aligned with the axis of the duct with 0 measured 
positive counterclockwise. The local system (£ 77 ,^) 

is Cartesian and is aligned with the reference (i.e., 
zeroth) airfoil. With no loss of generality, we assume 
that the 0 = 0 line and £ axes are coincident and 
along the leading edge of the reference airfoil (see 
Figure I). The relationship between the two systems is 
given by 



Figure 1 Cascade geometry and coordinates systems 

£ = jccos a s +r0sma, 

r\= -xsina, +r0cosa, (1) 


where a s is the stagger angle of the cascade. Eq. (1) 
holds as long as 0 is small (i.e., near the reference 
airfoil). The axial and radial coordinates are 
normalized by the tip radius, as arc the local 
coordinates. 0 = r H / r T is the hub-to-tip radius ratio 

with r H and r T denoting the hub and tip radii. In 
what follows, it is most convenient to describe the duct 
acoustic field in terms of the global coordinate system, 
while the building block of the cascade unsteady 
response is most easily derived in the local coordinate 
system. 

The representation of the duct acoustic field is 
outlined first followed by the construction of the 
cascade unsteady response model In both instances, 
we will take advantage of asymptotic methods (0 
simplif y the analysis. The connection between the two 
elemei ts is then established using the Green s function 
method. 
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Annular Duct Mode Expansion 


leads to separate equations in x and r given by 


For the duct mode acoustic field calculations, 
we assume the base flow to be axial. The propagation 
of acoustic pressure waves in such a system is 
governed by the convected wave equation which, in the 
duct coordinates, takes the form 


dr 


D 0 In' 

Or f 

= 0 

(2a) 

= iL + Iif 

d ) Id 2 

r l 1 

(2b) 

dx 2 r dr \ 

dr ) r 2 dO 2 

Si 

+ 

It 


(2c) 

t 

= 0 


(2d) 


I Duct 


p'{x y r,0J) is the acoustic pressure, a 0 the nominal 
speed of sound of the medium and V the base flow. 

Using the normal mode expansion approach, 
the most general form of p' with harmonic time 
dependence is given by 


CO 

p'(x,r,Oj) = Pmix.r) e ,{m0 ~ m,) (3) 

m- co 


p 2 — + 2ikM — + (k 


dx 


dx 


V) 


/»=o 


( 6 ) 


dr \_d^ 

dr 2 r dr 


+ /c 



(7a) 



= 0 


(7b) 


The term k 2 in these equations is the separation 
constant. We solve for the second equation first and 
postpone the solution of the first one until later in the 
analysis. 


Since Eq. (7a) is the Bessel’s Equation, it is 
usually solved in terms of the Bessel Functions of the 
First and Second kind of index m . The corresponding 
boundary value problem (i.e., Eq. (7b)) then leads to a 
transcendental equation for the eigenvalues of the 
annular duct (i.e., the permissible values of k ). The 
eigenvalues are generally obtained numerically by 
solving the transcendental equation iteratively. But 
here we pursue a different approach for computing 

K ’S. 

I r dr 

The integrating factor e ^ r = r 1/2 for Eq 
(7a) suggests the following change of dependent 

variable g m (r) = r~ vl g m (r ) . Rewriting Eqs. (7a & 
b) in terms of the new dependent variable yields 


where (o is the frequency and m the spinning (i.e., 
circumferential) mode order. Depending on whether 
one is concerned with tone noise or broadband noise, 
there will 'be a sum or an integral over all frequencies 
of interest. The form of pressure given by Eq. (3), 
simplifies Eq. (2a) to 




Krr l±) 

ax 2 rdr{ dr) 


nt~ d f2 

— r- + 2i/cM — + Ac 
r dx 



= 0 


(4a) 




(8a) 



(8b) 


< d_ 
^dr 



“r=<r.l 


= o 


(8c) 


k = ^ M — — , p= Vl-M 2 (4b) 

a () a o 

Here Ac is the reduced frequency and M is the Mach 
number of the mean flow . Introducing 

P' m (x.r) = f m (x)g m (r) (5) 


The distribution of the eigenvalues of the 
Bessel’s Equation is such that the smallest eigenvalue 
in the sequence is always greater than the index (i.e., 
Min(x) ^ m ) with the other elements of the sequence 
growing rapidly larger than m Therefore, as long as 
a> 0 (i.e., r H *0), w(r) will vary slowly with r for 

given k and m . This suggests that the WKB method 
(Olver [ 1974|) may provide a reasonable 
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approximation to the solution of Eq. (8a) and its 
eigenvalues. The accuracy of the approximation will, 
of course, depend on the behavior of »v(r), but as it 

will be shown the proposed approximation is 
surprisingly accurate, especially, for higher order 
eigenvalues as long as cr is not too close to zero. 


The asymptotic formula for the eigenvalues is 
derived in Appendix 1 The final result is 


nna 


tc— = 


+ J(ftxcr) 2 + <t(1 - a) 2 (2m 2 + 3 / 2 ) 
2cr(l - a) ~ 


(9) 



/? 


where k ^ is the non-dimensional eigenvalue and 
n = 0, l, 2, * • • its “radial” index set corresponding to 
each circumferential index m . The sy mbol = 
indicates that is an approximation to the exact 
eigenvalue. 

Eq. (9) is a simple formula but a remarkably 
accurate one as is shown in Figure 2 Here the exact 
(open bars) and approximate (solid bars) radial 
eigenvalues are plotted for two circumferential mode 
orders m - 6 and m - 30 for a hub-to-tip radius ratio 
of 0.5. The first 20 eigenvalues are shown in each 
case. Maximum approximation error is around 15% 
for low values of the radial index but as n increases 
the error diminishes rapidly. Given that the exact 
eigenvalues* must be determined numerically for each 
combination of m , n and o , the advantage of having 
an explicit formula for computing the duct eigenvalues 
accurately becomes apparent. 



n 


Figure 2 . Comparisons of exact and approximate radial 
eigenvalue distributions for two circumferential 
mode orders. For this calculation a = 0.5 . 

The prime in Eq. (10c) denotes differentiation with 
respect to the radial coordinate r and the sy mbol 91 
indicates the real part. 


The corresponding eigenfunctions, also 
derived in Appendix 1, are given by the expression 




i 


k r 


1- 


l-4m- 


x 91 


( . cr (*'mw- q ') 2ir.(g—.g) -iy.(r— .r) 


y a (K,r) = Kr + 


4m" -1 


8 Kr 

gj(k\ a) = 2a y/' a (*\ a) + 1 


(10a) 

(10b) 

(10c) 


The exact eigenvalues are solutions of the transcendental 
equation J'^K^a) Y' m ( k „„ ) - J' m ( k m „ ) = 0 . 


Eq. (10a) provides an efficient alternative to 
the exact * annular duct eigenfunctions which generally 
require recursive computation of the Bessel functions 
of various order. In Figure 3, a comparison of the 
exact and approximate eigenfunctions is shown. The 
first four normalized radial eigenfunctions are plotted 
for the circumferential mode m = 6 at a = 0.5 . 
Except for the lowest order approximate eigenfunction, 
w hich shows marked difference from the exact one, the 
rest arc consistently close to the exact eigenfunctions 
with the agreement improving as n increases. The 
fourth one, in particular, is virtually indistinguishable 
from the exact eigenfunction These results, and those 
shown n Figure 2, demonstrate that the formulas given 
by Eqs. (9) and (10a) provide useful approximation that 


The exact radial eigenfunctions are given by the expression 
= y (*• r) 

rrtft\ f m V mn ) , . . / v m \ mn / 
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(12a) 




Figure 3. Comparisons of exact and approximate radial 
eigenfunctions (normalized) for the circumferential 
mode order m = 6 with a - 05 Eigenfunctions 
corresponding to the first four radial orders are 
shown. 

can be used in place of the exact eigenvalues and 
eigenfunctions. 


f m (x ) > = 


,± , M , 1 2 

Jc ~ —k. + — r 

"" p 1 ' 


0 2b) 


where the plus sign indicates downstream running 
acoustic waves and the minus sign upstream running 
ones Combining the expressions for f m and g m , the 
general representation for the acoustic pressure field 
inside the duct is given by 


m co n l) 

03) 

For a given geometry' and operating condition, 
only a finite number of terms in the double infinite 
sums will constitute propagating waves. Waves (i.e., 

modes), for which k ^ is complex, are cut-off and do 
not contribute to the acoustic field. 

Annular Cascade Unsteady Response Model 

The complications introduced by the annular 
geometry of the cascade generally require numerical 
solution of the equations governing the unsteady 
pressure distribution on the cascade. To avoid a 
numerical treatment, the unsteady response will be 
constructed in a manner similar to that developed by 
Envia and Kerschen (1986) for the radiated noise from 
a rectilinear cascade. 


Eqs. (9) and (10a) together constitute the first 
of two main results of this paper. A detailed discussion 
of the behavior and analytic structure of these 
approximate eigenvalues and eigenfunctions is 
postponed to a future time 

With the solution of the radial eigenvalue 
problem determined, the most general solution of Eq. 
(7a) is given by 

oo 

= ( 11 ) 

n=0 

where C ^ s are acoustic mode amplitudes that will be 
determined later. 

Now, the solutions of Eq. (6) can be readily 
written down as 


Taking advantage of the weak coupling 
between the adjacent airfoils at high frequencies, these 
authors envisaged the cascade noise field as a sum over 
the individual airfoil fields with the blade-to-blade 
periodicity enforced. Here we modify the procedure to 
account for the annular geometry and deduce the 
unsteady pressure distribution on the cascade from the 
local acoustic field. The solution is of course 
approximate, but as was shown in the above reference, 
the approximation provides quite reasonable accuracy 
at high frequencies. 

The starting point for construction of the 
cascade response is the single airfoil solution. This 
solution is derived for the reference (i.e., zeroth) airfoil 
using the Wiener-Hopf technique in terms of “leading 
edge" and “trailing edge 11 responses that are combined 
to produce the response of the airfoil to incident flow 
perturbations. The details of the derivation, as well as 
a discussion of some sample calculations, are given in 
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Appendix 2 Here we present the final results 
expressed in die duct coordinate system. The solution 
is given by 


co (» CO 


= A u« cosf In - — — 


fj--ao y_0 /=0 


V \ — <T 


[G l (x,r.6)e‘ t ' [j ‘'* ) + G 2 (x.r,ey ri{ *'* i 1 (14a) 


r i (x,r.0)= zRj(x.r.O) 




l l \ 

~Y cosa, +rdsh\a s j (14b) 


Rj = ^x] +(r0)\ xj=x-(j-\)e 


(14c) 


where c is the airfoil chord The definitions for all the 
variables appearing in these equations can be found in 
Appendix 2. The triple summation provides the 
response over all possible gust modes that might exist 
inside the annulus as discussed in the beginning of 
Appendix 2. It should be emphasized that, at this point 
in the analysis, p' 0 is the unsteady pressure field in the 
vicinity of the zeroth airfoil and not just on its surface 
alone. After die cascade solution is constructed, the 
final expressions will be evaluated on die surface of the 
cascade. 


must hold between the responses of other airfoils and 
the response of the reference airfoil. That is to say, 

p' q {x,r.0.t) = p' 0 (x.r,O q ,t) e“ ,M ’ (17) 

Adding up the contribution of all airfoils we find 


p' c (x,r,0.t) = Y p' o (x.r,0 q ,t ) (18) 

q~ -oo 

where p' c is cascade response. The limits of the 
summation are set to infinity to ensure that the 
resulting solution and all of its derivatives are periodic 
around die annulus. Had we limited the sum to only 
N terms (corresponding to N airfoils), only the 
periodicity of the solution itself would have been 
guaranteed. Furthermore, as it will shown later, the 
infinite sum in Eq. (18) ensures that all circumferential 
modes that can exist in the duct are accommodated 

Taking advantage of the properties of the 
Dirac delta function, Eq. (18) can be rewritten as 





Pu(x,r,®j) e“ ,/Jf, S(e-0 q )de 

(19) 


Since the circumferential positions of any 
airfoil in the cascade is related to the circumferential 
position of the reference airfoil tlirough 

O q =0-q<p , q = 0, ±1, ±2, (15) 

where <p = 2n IN (A is the number of airfoils in the 
cascade), the acoustic response of any airfoil in the 
cascade can be related to the response of the reference 
airfoil. The required relationship can be established by 
noting that, on a per mode basis, the form of die 
incident gust on the q th airfoil expressed in die zeroth 
airfoil coordinates, is related to the form of the gust on 
the zeroth airfoil expressed in the^ th airfoil coordinate 
system through 


u, (x.r.Oj) = u _ (x 

'^airfoil q 0 th 




0™ airfoil 


(16) 


where we have taken advantage of the circumferential 
dependence of the gust as given by Eq. (A2.1) in 
Appendix 2. This suggests that an identical connection 


Interchanging the order of summation and integration 
and using the Poisson Sum formula (see Carrier et. al 
[1966] the sum in Eq. (19) can be converted to an 
equiva ent sum (say. with index t ) having die form 


p' c {x,r.0j) = 


_!_Y 

V2 ~7Up J 




* I Poix r^t) e ~ i(fi ' 2Ul9) ®d® (20) 


Given the complicated form of 

(see Eq. (14a)), the integral in Eq. (20) cannot be 
evaluated exaedy, but it can be computed 
asymptotically using the method of steepest descent in 
the high frequency limit (i.e., for k » 1 ). Substituting 
for p’ 0 \ jc.r,©,/) and grouping terms together, Eq (20) 
can be rewritten in canonical form as 
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where the phase functions contain alJ the phase 

terms that depend in the variable 0 and the amplitude 
terms A } include everything else in Eq. (14a). The 

phase functions are given by 




r0sin a, -{n + ltn / <p)S 


( 22 ) 


with Xj = x-(j - \)c The saddle points 0^ (i.e., 

the roots of (tVj I d<r) = = 0 ) are given by 


where 4^' is the second derivative of the phase 
function given in Eq. (22). 

The waves implied by Eq. (24a) are cut-off 
whenever the square roots in the expressions for 

^(© 7ip j (i.e., Eq. (24b)) become imaginary. In 

particular, consider the case for tone noise generated by 
the interaction of mean wakes of a 2?-bladed rotor 
with a V -vaned stator. For this situation the 
circumferential structure of the gust implies that 
H = pB where p is an integer If the stagger angle 

a s were zero, the cut-off condition for the lowest 
order spanwise mode (i.e., / = 0 ) would imply 



Carrying out the remaining steps in the 
analysis leads to the asymptotic approximation of the 
integral /, . The final result is given by 



(24b) 


pBOrr + (25) 

«oP r 

where we have taken advantage of the following facts: 
x\i =l) = k x M I f) = (or T I a n f ) , 
a) = pBQ . , 
q> — 2x IV 

and Q is the fan rotational speed. 

If the cut-off condition is satisfied for the tip 
radius, it will be satisfied for all other radii. Setting 
r = 1 in Eq. (25) and re-arranging the resulting 
equation, we arrive at the familiar Tyler-Sofnn cut-off 
criterion (see Tyler and Sofrin [1962)); 

>*L (26) 

pB fi 

where M T = Q/y is the tip Mach number of the rotor. 
Modes for which the above inequality holds are cut-off 
and do not contribute to the cascade response field 
For higher order spanwise modes and/or non-zero 
stagger angles the cut-off condition is somewhat more 
complicated. 
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Substituting for l t in Eq. (2 la), the final 
expression for unsteady pressure field in the vicinity of 
the cascade is given by 


CO OO CO to 


/J- -ao v~ 0 / 0 t- -CO 

\-{r Z Ki l CO * a s j Pt ) T ! 4 ^l( e i v )] 


o,(©, „)« 


JtK) 

g (© ■"«)* + T i( e ^)] 




iUn+2(x/e>)&-a/t] 

X e A 1 COS 


r - a 


In 
V 1 - <7 


(27) 


The last step in obtaining the unsteady pressure 
distribution on the cascade involves evaluating p' c on 
the airfoil surface. 


Eq. (27) is the second principal result of this 
paper. It provides a closed form expression for 
computing the asymptotic behavior of the unsteady 
pressure distribution on an annular cascade in the high 
frequency limit. The next step in the analysis involves 
relating the cascade pressure distribution to the duct 
acoustic modes. 


then converted to an equivalent description in the local 
airfoil coordinate system as given by Eqs (A2.2a 
through A2.2c). The resulting modal expansion will 
then be used as input to the cascade unsteady pressure 
formula (i.e., Eq. (27)). The result will be in turn 
integrated according to Eq. (28) to provide the modal 
pressure amplitudes C TO ’s. These will then be 
inserted in the expansion given by Eq. (13) to provide a 
description of the acoustic field inside the duct setup by 
the interaction of the incident gust with the annular 
cascade. 


Eq. (13) gives a description of the acoustic 
field in terms of duct modal pressure. However, such 
quantities as acoustic power can be easily calculated 
using (lie applicable formulas from Eq (13). 

In principle, the steps outlined above will 
apply whether one is interested in tone noise or in 
broadband noise. But, in practice, some reformulation 
of the formulas will be required for broadband noise 
calculations. In particular, since the description of the 
turbulent gust can only be given in terms of spatial 
correlations, the formulas derived in this paper must be 
recast in terms of expected values of acoustic power 
instead of acoustic pressure. The methodology, 
however, will remain exactly the same In an 
upcoming paper, detailed computations using the 
formulas derived here will be presented. 


CONCLUSIONS 


The methodology of Green’s function allows 
for the desired connection to be established between 
p' c and the values of (in Eq. (13)). The 

mathematical statement of this connection is given by 
(see, for example, Meyer and Envia [ 1996J) 



s 


where e is the unit normal perpendicular to the surface 
of the reference airfoil. The integration is carried out 
over the surface of the reference airfoil. Once C^, ’s 
are computed for each propagating inode according to 
Eq. (28), the duct acoustic field given by Eq. (13) is 
completely specified. 

In a typical computation, a description of the 
incident gust field inside the annulus will be developed 
first. This will be done according to Eq. (A2 1) and 


In summary, closed form asymptotic 
expressions for computing the noise field of an annular 
cascade inside a duct containing a uniform flow have 
been developed. By tailoring the analysis to the high 
frequency regime, numerical treatment of the equations 
involved has been avoided altogether 

There are two principal new results in this 
paper. The first is an alternative description of the duct 
acoustic modes that does not rely on the often-used 
Bessel function expansion. This results in simpler 
expressions for both the radial eigenvalues and 
eigenfi notions of the duct. In particular, the new 
representation provides an explicit approximate 
formula for the radial eigenvalues obviating the need 
for sclutions of the transcendental annular duct 
eigenvalue equation. Also, the radial eigenfunctions 
are represented in terms of exponentials eliminating die 
numerical problems associated with generating the 
Bessel functions on a computer 

The second new element is the construction of 
an unsteady response model for an annular cascade 
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4 m‘ - 1 


using the single airfoil solution as the building block. 
The new construction allows for the boundary 
conditions on both the cascade and duct walls to be 
satisfied simultaneously adding a new level of realism 
to the noise calculations. Preliminary results were 
presented which demonstrate the utility and 
effectiveness of the new elements 

Cascade discrete tone and broadband noise 
predictions using the new F formulas will be presented in 
a future paper along with comparisons with results 
from simpler 2D-strip based models and rectilinear 
cascade models. 


APPENDIX 1 






8 a: r 


(A 1 .3) 


where the subscript a indicates that Eq. (A 1.3) is an 
approximation to Eq. (A 1.2). The constant term in the 

expansion (i.e., the ^( x ‘°) term) is ignored since it 

only contributes a multiplicative constant to the 
expression for g m and therefore does not alter the 
subsequent results. 

Substituting y/ tt for y/ in Eq (Alla) and 
noting that [**>(/■)] '^ = K~' n +o(/r" 2 ) , the boundary 
condition at r = o determines C 2 in terms of C } via 


Asymptotic Solution of the 
Radial Eigenvalue Problem 

The WKB approximate solution of Eq. (8a) is 

given by 


C 2 = C, 


gr(*r,cr) ^ 2 iy m (K,a) 
gt*(a\ cr) 


(A1.4a) 


ot(a",cf) = 2<7^(a',<t) + i (A 1.4b) 


S m (r) = 


^lMr) 


i K'(r) dr -i I 1 

C } e j { +C 2 e J 


’(') dr 


(Alla) 


subject to 


where the prime denotes the derivative with respect to 
r and the asterisk a complex conjugate quantity 
Enforcing the boundary condition at r = 1 and 
substituting for C 2 from Eqs. (A 1.4a & b) yields 



(A 1.1b) 


1) *«(*•**)] _ QT (at, l) ^ (AT. Cj) 

xb{k, l) cr) 


(AI.5) 


where C } and C 2 are arbitrary constants and the 
symbol = denotes the approximate nature of the 
solution. The integration indicated in the phase of the 
exponential terms can be carried out explicitly, 
yielding 


y'(A' s r) = ^j > w(r) dr = -^4(a : r)“ - 4 + 1 

w 

1 + +1 

V 1 - 4m 2 


- — yj \ - 4 w { 

2 


- log^ yj\-4nt 2 + 2 log(2) 


(A1.2) 


Assuming that k » 1 , Eq. (A1 .2) can be simplified by 
expanding it at that limit. The simplified phase 

function to o[k ~ 2 ) is given by 


which is the equation that determines the non- 
dimensional eigenvalue k . In its present form, Eq 
(A 1.5) is transcendental in k and not amenable to an 
exact analytical solution. But if we take the logarithm 
of both sides and expand the resulting equation to 

consistent with the level of approximation 
introduced in obtaining y / a , we find 

2[^ fl (AT, \)-y/ a {K. cr)] = ^——^+2nn (A1.6) 

where #i = 0, 1, 2. ••*. In view of Eq. (A1.3), Eq. 
(A 1.6) is quadratic in k and can be readily solved to 
yield 


/i/r<7 + 


= ■ 


^(/f/rcr) 4 ' + a (\ - cr)^2m 2 + 3/2) 
2cr(l- cr) 


(A1.7) 
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where only the positive root has been retained since 
only the positive eigenvalues are relevant here. Note 
that, the pair of indices m.n denotes a doubly infinite 
set. Eq. (A 1.7) is a remarkably simple formula for the 
approximate eigenvalues of the annulus. The accuracy 
of this formula is investigated in the analysis section of 
this paper. 

Inserting the expressions for C 2 and k m in 
Eq. (Al la), we find the most general representation of 
?«('■) as 


gm(r) = 


r 

u 


1-4 m~ 


x 91 


-..»•) j tg ( y m«. <T ) g 2»y.| 




(A 1.8) 


where denotes the real part of the quantity and 
is the constant C \ renamed. To improve the accuracy 
of the eigenfunction approximation the next order term 

r i *1/2 

in the expansion of [^(r)J “ has also been kept 

resulting in the term in front of C ^ 


APPENDIX 2 


where ’s are the radial eigenfunctions of the duct 

(see Eq. (10a)). The axial wavenumber is given by 
k x -k! M owing to the requirement that the gust be 
convected by the mean flow (Recall that k = (or T fa 0 ). 
A pv 's denote the gust modal amplitudes. 

Along the leading edge of the reference airfoil 
0 is zero. Setting 0-0 in Eq. (A2. 1), expressing the 
resulting equation in the local coordinate system, and 
expanding the £ dependence in a Cosine Fourier 
series yields 


00 OO OO 


*r»r? 


//- - oo k-0 / - 0 

x cos(lnC) (A2.2a) 

r 4 = k x cosa s , y T} = -k x sina s (A2.2b) 


i 


a ^= A „v I <*v[cr + (l-a)C]cos(/<)rfC 


(A2.2c) 


Eq. (A2.2a) provides the general modal representation 
of the incident gust in the local coordinate system 


Asymptotic Solution Of 

The Single Airfoil Unsteady Pressure Distribution 

It is mathematically more convenient to solve 
this problem in terms of the acoustic velocity potential 
(denoted here by $ ) instead of the pressure. Once the 
solution is obtained, the pressure can be calculated 
from ^ via the linearized momentum equation, i.e., 

p* = vv ^ ere P 1S the ambient density. 

As was indicated earlier, the unsteady 
pressure response will be formulated and solved in the 
local coordinate system (£ It is, therefore, 

necessary to represent the incident gust in these 
coordinates. The most general representation, in the 
global coordinate system, of the relevant component of 
the convected gust is given by 

OD 00 

fj - -- oo y=0 

(A2.1) 


In view of Eqs. (A2.2a) through (A2.2c), the 
mathematical statement of the unsteady response of a 
single airfoil subject to a single incident gust mode is 
given by 

[v 2 -flo 2 -^W.^C,0 = ° (A2.3a) 

^ =-e { ' rt *~ m, \os (lx£) on 0<4<c (A2.3b) 

dn 

9 

/>'£“ = 0 on (A2.3c) 

where c is the airfoil chord The transverse wave 
number y n does not enter Utese equations owing to the 

form cf the boundary condition in Eq. (A2.3b). The 

...... , D 0 d , a 

convective derivative is given by — 1 - = — + K, — 

Dt dt 5 

where V^-U cos a, is the speed of the uniform base 
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flow, which now is aligned with the airfoil chord. The 
last equation enforces the continuity of the acoustic 
pressure downstream of the trailing edge. These 
equations are supplemented by the Sommerfeld 
radiation condition at infinity and the continuity’ of 
derivatives of $ off the airfoil. 


A Fourier transform in reduces Eqs. 
(A2.6a) to an ODE whose general solution is given by 

=sgn(» 7 ) C,(A)<? ^ x ^ (A2.7a) 


The exact solution of Eq. (A2.3a), subject to 
all of the specified boundary conditions, can not be 
obtained in closed form But, it has long been 

established that the solution may be developed as a 
convergent series ^ + <0 2 + " ,n which each term 
is a solution of Eq. (A2.3a) but with one of the two 
edges moved off to infinity. As a result, for each term 
in this sequence, some of the boundary conditions can 
be ignored leading to a simpler problem which can be 
solved in closed form. Landahl [ 1958J has shown that 
the series converges for all frequencies with the rate of 
convergence increasing rapidly as the frequency is 
increased. It will be shown later that, for our purposes, 
the first two terms in the sequence provide sufficient 
accuracy. 

The Leadine Edge Solution 


d+ 1 

dr] 


rf=0 


^f2x 


Pi 


(A2.7b) 


where X is the Fourier transform variable, the 


transform of and Eq. (A2.7b) the transform of Eq. 
(A2.6b). The branch cuts of the square root in Eq. 
(A2.7a) are chosen so that the real part of the root is 
always positive. The unknown constant C,(A) can be 
expressed in terms of known functions by substituting 
for in Eq. (A2.7b) and taking advantage of the 
requirement that the acoustic velocity be continuous 
everywhere. The resulting equation can then be easily 
solved through the use of the Wiener-Hopf technique 
(see Noble [1958]). The final expression is given by 


The leading edge solution must satisfy Eq 
(A2.3a) subject to 


= -e^ r< ^ "^cos(fjHj') on 0<£<ao (A2.4) 

dT > „ 

The introduction of the following change of variables 



(A2.5a) 

Pz a o 

(A2.5b) 

simplifies Eqs. (A2.3a) and (A2.4) to 


( d 2 d 2 2 )- n 

(A2.6a) 

W 2 drj 2 J 



c, (*) = -- 


Vi 


(A2.8) 


+ X 


X + 


n 

Pi) 


yft + Z 


Substituting for C,(A) in Eq. (A2.7a) and applying an 
inverse Fourier transform to the resulting expression 
yields the solution for . Reverting to the original 
dependent variable and computing the acoustic 

pressure (i.e., p[ = ) 8 ives 




-ipy S8n(l),. cos(<<) 


x e 




fi+x 


dX (A2.9) 


Qn o < < oo (A2.6b) 

dl > 



(A2.6c) 


The pole in Eq. (A2.8) does not appear in the 
expression for pressure since it represents a solenoidal 
(i.e , pressure-free) velocity. The integral in Eq (A2.9) 
can be evaluated explicitly in terms of elementary 
functions. The development is straightforward and, 
therefore, only the final answer will be given here The 
solution is 
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Pl(4‘.v^t) = - 


p V sgn(i;) e 


ijr/4 


. -COS 

1 *■ -*cos (IxC) (A2.10a) 


Following the same procedure that was 
outlined for the leading edge solution, first a reduced 
potential is introduced, 

* 2 =* 2 el J cos (/<) (A2.13) 


u' = tan 


\ JL 


(A2. 10b) 


and then a Fourier transform in (i.e., £' = £//?,*) is 


Eq. (A2.I0a) gives the pressure field 
associated with the leading edge everywhere. Note that 
the pressure has a square root singularity at the leading 
edge and decays algebraically away from it. The 
airfoil chord, c , does not appear in the leading edge 
solution since the trailing edge was moved off to the 
downstream infinity. The pressure distribution on the 
airfoil (i.e., for u’ = 0 ) is given by 


O.C0 


_ p V e'"' 4 e 

= +■ 






-cos (/<) (A2.ll) 


The plus sign indicates the pressure on Hie upper 
surface and the minus sign the pressure on the lower 
surface. Since the leading edge solution cannot satisfy- 
the Kutta condition at the trailing edge (i.e., p{ is not 
continuous there), the trailing-edge solution must be 
added to p\ to rectify this shortcoming. 


applied allowing the formal solution for (the 
Fourier transform of ^ 2 ) to be found. The result is 

=sgn(> 7 ) C 2 (X)e v^W (A2.14a) 


subject to the transformed version of Eq. (A2. 12b), 


Ps) 


n=o* 


= 0 on £>0 (A2.14b) 


rf=0 


where the terms on the left-hand side are the Fourier 
transforms of the acoustic pressures associated with 
and respectively, and c'= c ! . 


The Trailing Edge Solution 

Now we let the leading edge move off to the 
upstream infinity and place the origin on (lie trailing 

edge. The new chordwise coordinate £ is simply 

£ = The trailing edge solution must also 

satisfy Eq (A2.3a), but the appropriate boundary 
conditions are now given by 


The integral in Eq. (A2.14b) can be evaluated 
explicitly in terms of the Fresnel Integrals. But, since 
we are only interested in the pressure distribution 
upstream of the trailing edge, it suffices to determine 
the behavior of the integral near the trailing edge (i.e., 

for £ # = 0 ). Expanding the denominator of the 
integraid and integrating term by term yields the 
power series 


d# 2 

Or) 


VI 

c 

o 

© 

II 

A2.12a) 1 

17=0 





. VI 


d? = 


*’+c’ 


'&{*■ + X) 2<[c^(X + xf 


7 _ 

= 0 on £>0 (All 2b) 

rr<> 

where p{ in Eq. (A2. 12b) is the leading edge solution 
(i.e., Eq. (A2.ll)) expressed in the shifted chordwise 




+0({A+x) i ) (A2.15) 


of wh ch only the leading term is kept in the 
subsequent analysis. This approximation simplifies the 
algebra significantly while still guaranteeing that the 
Kutta condition is satisfied at the trailing edge 


coordinate £ . 
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In a slightly more complicated Wiener-Hopf 
analysis, the unknown constant C 2 (> l) in Eq. (A2.14a) 
is found to be 


C 2 (A) - -- 




n4c* , 



( \ 


Pi) 


(*+z) -Mz 

(A2.16) 


representation of the unsteady pressure distribution on 
the airfoil. The result is 

ri(f, o.c.1). 

_ e ‘ll* (>•<*<■".))«•■ -'l cos((<) 

<A2m 


With C 2 (/l) known, Eq. (A2.14a) can be formally 
inverted. Switching back to the original dependent 
variable and computing the acoustic pressure we 
find 




i[X 


pA 2? 


+ Z 


xe -N*.^.)i--'] C0S ;(l<) 


i: 


" p 

if*t 

e -iit'-ypl 2 -Z 2 \r?l 

yj x ~Z ( X + x) 

(MM) 


dl 


The above integral can be evaluated explicitly. After a 
fair amount of algebra, the final result, expressed in the 
leading edge coordinates (£', 7 ,^) , is given by 


= P ^sgn(>y) e ixl4 CQS ^^ 

J*Pf c ' Jrs+P 4 z 


2e 


-ixt 4 


^' + ^ 




s 4 .= 


0 if £'-c '<0 

1 if -c* >0 


(A2. 18a) 
(A2. 18b) 


where 


fM-JV 


r/a is the Fresnel Integral. Note 


that, unlike the leading edge solution, the trailing edge 
solution is not singular 


Evaluating p\ on the airfoil (i.e., 7 = 0 ) and 
combining it with Eq. (A2.ll), yields the two-term 


The combined solution satisfies the Kutta condition at 
the trailing edge (i.e., at = c * ) since 

F( 0 ) = 4ne ix ' A 1 2 . 

In Figure 4 magnitude of the unsteady surface 
pressure as predicted by Eq. (A2.19) is plotted for 
typical (high) reduced frequency and Mach number 
values. The corresponding plots for Pi' ± (£\0,£\f) and 

P 2 t (£\0,£,f) are also shown for comparison. For this 

calculation the gust is assumed to be uniform along the 
span. Therefore, Figure 4 effectively shows the 
chordwise pressure distribution. 



Figure 4 Typical chordwise unsteady pressure distribution 
on a single airfoil. p\ and p\ are leading edge 

and trailing edge solutions, and p ' the combined 
solution (x£,^) = (lO,0) and = 0.5 for 
this calculation. 
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Figure 5. Predicted real and imaginary parts of the 3D 
unsteady pressure distribution on a single airfoil for 
a gust with radial structure The gust has the form 

Ug - <I> 62 (r) e 7 ^ » with Y t = 22 and M ^ - 0.5 

Note that p\ ’s primary contribution is in the 
vicinity of the trailing edge where it enforces the Kutta 
condition by canceling out the leading edge 
contribution there. Everywhere else it is only a small 
correction to the leading edge solution. The third term 
in LandahTs series will be correspondingly smaller in 
comparison with the leading edge solution. At high 
frequencies, therefore, the additional complexity 
incurred by adding the third term is hardly justified for 
a very' small improvement in the accuracy of the 


Re(p') 



Figure 6. Predicted real and imaginary parts of the unsteady 
pressure distribution for a single airfoil using 21)- 
strip approximation The gust parameters are the 
same as those in Figure 5 

solution. Hence, we shall truncate the series at this 
point and use the two-term representation in the 
subsequent analysis. 

In Figure 5 the predicted real and imaginary 
parts of the unsteady pressure distribution for a single 
airfoil as given by Eq. (A2.19)) arc shown for a more 

realisti; gust distribution given by u g -<\y 62 {r) e r ^ 
with y ^ - 22 , a - 0.5 and M ^ = 0.5 . The radial 
structure of the gust means that there will be a 
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spectrum of spanwise wavenumbers as dictated by Eq. 
(A2.2a) and their associated spanwise amplitudes as 
given by Eq. (A2.2c). Only the first five terms in that 
series were taken into account for this calculation since 
the higher order terms result in cut-off waves (i.e.,^ . 
given by Eq. (A2.6c), is complex for / > 5 ). The 
corresponding amplitudes of the spanwise modes are 
given by (0.067, 0.127, 0.375, -0.028, -0.004). For 

the sake of comparison the corresponding 2D-strip 
approximation to the unsteady pressure distribution is 
plotted in Figure 6. For this calculation, the variation 
of the gust amplitude along the span was taken into 
account but not its spanwise wavenumber spectrum 
While overall structure of the 3D and 2D-strip 
distributions look similar, they exhibit differences in 
their detail. A discussion of the importance of these 
differences and their effect on the noise fields is 
postponed to a later paper 

In order to utilize the formulas developed in 
this appendix in the analysis outlined in the report, it is 
necessary to express them in the duct coordinates. 
After some simplification and grouping of terms, the 
final expressions are given by 


oo <r> oo 


p’ ltj (x.r,0j) = ^ ^ A M «Gj{x,r,0)co^lxj-^ 

fj- -<n O t ~ () 


x e 


r .( x.r.0) oft 


(A2.20a) 


T,(x-r0)= zRj(x,r,0) 


Mi 


-y ^ — — (jccosa, +rOsina s ) (A2.20b) 

Pi 


,9= tan 1 


-jcsina, +r0cosa, > 
jrcosa, +r0sina f y 


(A.20g) 


where j = 1 gives the leading edge solution and j = 2 
the trailing edge solution. Eq. (A2.20a) represent the 
unsteady pressure field in the vicinity of a single airfoil 
bound by the walls of the duct with the index " 0 ” on 
p' 0j indicating that the solution applies to the zeroth 

airfoil. 
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